Quantum Zeno-based control mechanism for molecular fragmentation 
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A quantum control mechanism is proposed for molecular fragmentation processes within a scenario 
grounded on the quantum Zeno effect. In particular, we focus on the van der Waals Ne-Br2 complex, 
which displays two competing dissociation channels via vibrational and electronic predissociation. 
Accordingly, realistic three dimensional wave packet simulations are carried out by using ab initio 
interaction potentials recently obtained to reproduce available experimental data. Two numerical 
models to simulate the repeated measurements are reported and analyzed. It is found that the 
otherwise fast vibrational predissociation is slowed down in favor of the slow electronic (double 
fragmentation) predissociation, which is enhanced by several orders of magnitude. Based on these 
theoretical predictions, some hints to experimentalists to confirm their validity are also proposed. 



The quantum Zeno effect (QZE) was theoretically pre- 
dicted in 1977 by Misra and Sudarshan 1 and experimen- 
tally observed by Itano et al£ almost 15 years later. 
The opposite phenomenon, the so-called anti-Zeno effect 
(AZE), was later on predicted^ and observed immedi- 
ately afterwards. 5 These are well-known phenomena in 
the literature of quantum theory of measurement and 
open quantum systems, mainly in connection with deco- 
herence processes^ If a series of successive observations 
or measurements are carried out on an unstable system 
at very short time intervals r, its decay can be either in- 
hibited (QZE) or accelerated (AZE). These two remark- 
able behaviors with time become apparent when looking 
at the survival probability Pit), namely the probability 
to find the system in the same initial state, \tp(0)), at a 
certain subsequent time t, 

p{t) = |<v(o)|^)>| 2 = I Wo)| e -^/^(o))| 2 . (i) 

Here, H = Ho + V encompasses the system free-evolution 
Hamiltonian, Hq, and the system interaction, V. Because 
of the continuous measurement (or projection) process, 
the system decay can be described, in general, as^ 



P(t) 



(2) 



where 7(t) is the measured-modified decay rate, which 
depends on t. This rate is different from the free- 
measurement system decay rate, 70. The decay rate 7(7") 
can be expressed^ in terms of the spectral density of final 
states, G(w), and the measured- induced initial state level 
width, F(u>; r), as 



(3) 



7(t) = 2tt / G(u)F{uj;T)duj. 



The measurement frequency v = 1/r is related to the ini- 
tial state energy uncertainty, AE, as AEjv ~ h. Hence, 
the decay rate is essentially determined by the spectral 
density profile within a bandwidth v around its energy 
level. According to Kofman and Kurizki,' 3 Eq. con- 
stitutes a universal result: frequent measurements on a 
given initial state generally lead to its dephasing through 
randomization of the corresponding phase. Depending 
on whether v is much larger or much shorter than the 



spectral density width with respect to the center of grav- 
ity of G(uj), either QZE or AZE is expected, respec- 
tively. In this regard, ID systems have been largely con- 
sidered theoretically^^ while multidimensional realistic 
systems are almost unexplored. Furthermore, Kofman 
and KurizkiiS. have already discussed the idea of using 
QZE to control the decay rate of quantum systems cou- 
pled to a continuum reservoir by means of weak pertur- 
bations. The experiment carried out by Fischer et al., 5 
with cold atoms initially trapped in an optical lattice, 
precisely shows the appearance of QZE and AZE by re- 
peatedly on/off switching of the coupling between an al- 
most bound state and a continuum. On the other hand, 
it has also been shownii that periodic coherent pulses 
acting between the decaying level and an auxiliary one, 
can either inhibit or accelerate the decay into a reser- 
voir. In chemical reactivity, evidences of AZE have been 
shown by Prezhdo, 12 who also considered this mecha- 
nism as a novel route to quantum controlling in chem- 
istry. This is a step towards the control of decoherence 
and entanglement in molecular systems which can handle 
large amounts of quantum information^ - — 

In this work, on the contrary, a mechanism is proposed 
to quantum controlling the time-evolution of several com- 
peting fragmentation processes of an isolated molecule in 
gas phase, based on the QZE. Thus, by doing repeated 
measurements on one of the dissociation channels the 
outcome of an unimolecular fragmentation can be con- 
trolled. Isolated molecular systems allow to interpret 
experimental results based on realistic multidimensional 
models, allowing accurate quantum dynamical studies. 

To this end, rare gas-halogen diatom van der Waals 
complexes (Rg— X2) are ideal systems displaying a nat- 
ural separation between the fast intramolecular X2 vi- 
bration and the slow intramolecular Rg— X2 motions. 7 
The van der Waals (vdW) complex is prepared in the B 
excited electronic and vibrational v state, Rg— X 2 (-B, v), 
by the optical B <~ X electronic transition. The vibra- 
tional energy, initially deposited in the X2 subunit, flows 
towards the weak Rg— X 2 mode producing the fragmen- 
tation through vibrational predissociation (VP) , 

Rg-X 2 (B,v)^Rg + X 2 (B,v' <v). (4) 

The dissociation process is monitored by detecting the 
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FIG. 1: Schematics of the VP/EP fragmentation process 
considered here, including the pump/probe and Zeno- type 
pulses (see text for details). In the inset, measured^ and 
simulated^ total lifetimes. 

X2 (£?,?/) fragments by a second probe laser. The pres- 
ence of the Rg atom induces electronic couplings in the 
X2 subunit. This may eventually produce electronic pre- 
dissociation (EP) via double fragmentation, 

Rg-X 2 (B,v) ^Rg + X(P3 /2 )+X(P 3/2 ), (5) 

indicated schematically in Fig. [TJ This process was first 
detected for Ar-l2 complexJ^r— The Ne-Br2 system con- 
stitutes a benchmark since the existing VP/EP compe- 
tition has been extensively studied. Numerical results 
obtained using 3D potential energy surfaces issued from 
ab initio calculations- 6 have been found to be in excel- 
lent agreement with experimental complex lifetimes^ as 
seen in the inset of Fig. [T] In this figure, oscillations are 
attributed to the EP process due to the Frank-Condon 
factors between the electronic B state and the dissocia- 
tive 2 g and C states of Br 2 . As expected from the energy 
gap law,— the VP rate increases monotonically with v; 
for v > 15, it becomes more efficient and dominates. 
Based on the excellent realistic model describing the ex- 
perimental results, here we study the NeBr2(-B,v =20) 
decay using 3D wave packet (WP) simulations under fre- 
quent measurements. As is shown below, this allows us to 
slow down the otherwise dominant VP for v=20 in favor 
of the slower EP, which involves two different dissociative 
electronic states, the 2 g and C states of Br2. 

Let us first only consider the VP process for Ne- 
Br 2 (B,u = 20) in order to illustrate the QZE and AZE. 
This is carried out by expanding the WP as 23 

\^(t)) = J2^(R,0,t)ip v/ (r). (6) 

v' 



Here, (r) denote the vibrational eigenfunctions of the 
Br 2 (£?) subunit; & v /(R,9,t) depends on the distance R 
between Ne and the Br 2 center-of-mass and the angle 
9 between the two Jacobi vectors, r and R. The pre- 
dissociation lifetimes are of the order of 10-100 picosec- 
onds, while the quadratic behavior of the decay occurs 
at a much shorter time scale (of the order of femtosec- 
onds or shorter). The inclusion of a new pulse (Zeno 
pulse in Fig. [T]) in the propagation is computationally 
very demanding, and a first order perturbation approach 
is adopted for simplicity to deal with the corresponding 
absorption. Within this approximation, the effect of the 
Zeno pulse can be described by two limiting cases. First, 
the population depletion method, where the pulse fully 
promotes the dissociative $ 1 ,_i(_R, 9, t) component of the 
WP to a higher electronic state by optical pumping. This 
is simulated by projecting the WP on the different vibra- 
tional states and then removing the v — 1 component at 
each measure time nr. Second, the dephasing or ran- 
domization method where the measurement introduces a 
random phase in the & v -i(R,9,t) component as a per- 
turbation, simulated by multiplying $ t ,_i(i?, 9, t) by a 
random phase depending on the coordinates R and 9. 

The effect of these two measurement schemes on P(t) 
can be explained by formally expanding the WP as 

|*(t)> =a b (t)e- iE ^ h \b) + J2af(t)e- iEft/h \f), (?) 

/ 

where \b) and |/) denote the initial bound state and fi- 
nal dissociative states of Ho, respectively. At very short 
times, the decay of \b) into some final |/) is described by 
first-order perturbation theory as^ 

da b i T/ -i(E b -E f )t/h io\ 

By means of a series of repeated measurements at time in- 
tervals r — > 0, it is possible to freeze the evolution of the 
initial state, slowing down the transition by destro ying 
the coherence on the right-hand side of Eq. ([8]) (Ref. \2at) 
(for each (R,9) pair, there is a given positive/negative 
sign for the probability amplitude and the integrated am- 
plitude vanishes). As seen in Fig. [2] the typical QZE 
and AZE behavior is numerically reproduced for Ne- 
Br 2 (-B, v=20) using the randomization model with differ- 
ent values for r. With the depopulation method nearly 
indistinguishable results are also obtained, demonstrat- 
ing the robustness of the approximations made. When 
considering repeated measurements with short r, the life- 
time increases by an order of magnitude, from 27 ps in 
the free (unperturbed) case to about 320 ps in the case of 
r=5 fs. A similar deceleration is also found for r < 50 fs. 
This is a clear manifestation of the QZE. For r = 100 fs, 
the lifetime is 20 ps, approximately the same as in the 
free case, while for r > 100 fs, VP is enhanced, yielding 
a lifetime of 13 fs, according to the AZE. The exponen- 
tial decays for different values of t fulfill the predicted 
behavior given by Eq. @. In Fig. [5] the fittings to that 
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FIG. 2: Survival probability as a function of time for the 
free (unperturbed) evolution (red) and using different time 
intervals r for the randomization procedure. In the inset, en- 
largement on the short-time dynamical regime. Open circles 
are fits to the exponential decay of Eq. [2] with 70 = 0.098 
(free) and 7 (t) = 0.008 (r=5 fs), 0.13 (r=100 fs) and 0.199 
(r=200 fs), all given in cm -1 . 



expression are shown with open circles. The QZE is only 
found for short values of t, where the decay is quadratic. 
The Zeno behavior is lost after 0.1-0.2 ps, just when the 
curvature of the free population changes. It should be 
noticed that the v — 1 dissociative component of the WP 
is measured instead of the initial v component. How- 
ever, the usual QZE and AZE behaviors are found. This 
is easily understood within the randomization measure- 
ment model. The dissociation in this case is slowed down 
by introducing a random phase difference between the 
bound and dissociative components of the WP. Such a 
difference can be applied either to the initial, v, or dis- 
sociative, v — 1, component, leading to the same results, 
when only one fragmentation channel is considered. 

In the case of several fragmentation channels, mea- 
suring the initial state would slow down the total decay 
towards all the dissociation channels. On the contrary, 
measuring one of the dissociative states, would slow down 
selectively a given dissociation channel leaving nearly un- 
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FIG. 3: Electronic populations of the 2 g and C dissocia- 
tive states of the Br2 as a function of time, due to EP for 
the free dynamics and under repeated measurements with 
t=5 fs. Solid lines correspond to numerical calculations. 
Open circles correspond to fittings of the form Qa(l-e _1 ' T ''), 
with 70=0. 106 cm" 1 , 7 (r=5 fs)=0.0083 cm" 1 . The branch- 
ing ratios are Q° 2a =0.0133, Q% =0.0074, Q% =0.325 and 
Qc =0.103. 

changed the other ones. This scenario is precisely the 
one used here when considering the competition of the 
VP channel described above together with the EP chan- 
nel towards the 2g and C dissociative states of Br2. It 
is therefore expected that the slow electronic fragmen- 
tation described by Eq. (JSJ can be enhanced by slowing 
down the dominant VP channel. To deal with the two 
dissociation mechanisms, the total WP is written as2£ 

\9(t)) = Y,*°(r,R,0,t) X a, (9) 

a 

where \ a (a — B. 2 g and C) is the electronic wave func- 
tion and & a (r, R, 6) is expanded on a grid of 256 x 160x40 
points for the three Jacobi coordinates.*™ Due to the large 
grid involved, the calculations have been parallelized in 
the angle 9^ The randomization model employed above 
becomes unstable numerically for long times. For this 
reason, the measurements are simulated by the depopula- 
tion method, by projecting the $ a =s(^) R, 0) component 
onto the initial state, removing all other possible con- 
tributions describing the VP dynamics at each t. The 
typical QZE oscillations are again observed as in the pre- 
vious case. The slowing down effect introduced by the 
repeated measurements is again remarkably reproduced 
and rather important. Since the VP process is slowed 
down, the whole lifetime increases. 

Under such conditions, EP rates become larger or of 
the same order than VP ones. Thus, the electronic pop- 
ulation for the two dissociative electronic channels in- 
creases. In Fig. |3] the electronic population of the Br2 
2 g and C dissociative states, is displayed as a function of 
time. An analytical fitting to the exponential law is also 
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shown; the electronic branching ratios are obtained from 
this fitting. The branching ratio of the 2 g state changes 
from 1.3% in the free case to 32.5% with repeated mea- 
surements every t=5 fs; for the C state, this ratio changes 
from 0.74% to 10.3%. This implies the branching ratios 
of these two states increase by factors 25 and 14, respec- 
tively. This is a clear evidence of how the Zeno pulse 
can slow down the otherwise dominant VP, producing 
an important enhancement of EP, which leads to a total 
Br+Br+Ne fragmentation. 

For the particular case studied here, experiments could 
be performed by simply adding a third pulse, the Zeno- 
pulsc, to the usual pump-probe experiment^ The fre- 
quency of this Zeno-pulsc should be tuned into the B — > 
E transition maximizing the absorption intensity. It is 
expected that the intensity should be high to promote 
the population of the (B, v — 1) channel as much as pos- 
sible. The Zeno-pulses should be short, of the order of 
1 fs in order to allow a measurement frequency v of the 
order of 0.2 fs -1 . The relative phase of the pulses are 
not expected to play a major role since the (coherent or 
incoherent) dynamics on the upper E state is not of in- 
terest here. These effects should hold for r — > 0, and it 
can be envisaged to use a continuous laser to generate 
the Zeno pulse, as recently proposed by Kurizky and co- 
workers 27 in a different context. In this way, the direct 
detection of EP products, namely the atomic Br( 2 P 3 / 2 ) 
fragments would allow to measure EP rates directly in 
absence or by diminishing VP. Thus, the oscillations of 
EP rates as a function of v associated with the electronic 



B/2g or B jC Frank- Condon factors could then be better 
interpreted. This could help to unravel the complex dy- 
namics in other similar systems, such as the prototypical 
Ar-l2i2£ This working scheme will also allow the other 
way around, that is, to enhance VP in cases where EP is 
dominant, such as for v=\\ in NcBr2. 

In summary, the QZE-based control mechanism is of 
rather general applicability, and only requires two com- 
peting decay channels. By acting on one of them, by 
slowing down or accelerating it using the QZE or AZE, 
respectively, one can obtain a change of the branching 
ratios of one or two orders of magnitude. This quantum 
control method is specially well adapted for slow decay 
dynamics, where other control methods requiring coher- 
ence fails, such as the coherence control method. When 
including the Zeno-pulse in the simulations, there may be 
an attenuation of the processes described, but the results 
obtained are expected to be rather robust and will not 
change the major physics here proposed. 
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